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Abstract 

We construct operators that describe power corrections in mixed collinear-ultrasoft processes in 
QCD. We treat the ultrasoft-collinear Lagrangian to O(A^), and heavy-to-light currents involving 
collinear quarks to 0{X) including new three body currents. A complete gauge invariant basis is 
derived which has a full reduction in Dirac structures and is valid for matching at any order in 
as- The full set of reparameterization invariance (RPI) constraints are included, and are found 
to restrict the number of parameters appearing in Wilson coefficients and rule out some classes 
of operators. The QCD ultrasoft-collinear Lagrangian has two O(A^) operators in its gauge in¬ 
variant form. For the 0{X) heavy-to-light currents there are (4,4,14,14,21) subleading (scalar, 
pseudo-scalar, vector, axial-vector, tensor) currents, where (1,1,4,4, 7) have coefficients that are 
not determined by RPI. In a frame where u_l = 0 and n-v = 1 the total number of currents reduces 
to (2,2,8,8,13), but the number of undetermined coefficients is the same. The role of these oper¬ 
ators and universality of jet functions in the factorization theorem for heavy-to-light form factors 
is discussed. 
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I. INTRODUCTION 


The soft-collinear effective theory (SCET) constructed in |l|, l^y, offers a systematic 
description of processes involving energetic particles. It has an expansion in a small param¬ 
eter A ~ p^/Q, where p± is a typical transverse momenta and Q the large energy scale. 
Hard exclusive and inclusive processes in QCD are usually described using the powerful 
techniques of QCD factorization and light-cone expansions jl,0- SCET encompasses and 
extends these frameworks, and in particular allows a model independent description of ef¬ 
fects caused by the interplay between energetic collinear particles and soft particles beyond 
leading order in the power expansion. These effects can be described in a rigorous way based 
solely on QCD, but are not included in purely collinear expansions. The study of opera¬ 
tors that describe these mixed collinear-ultrasoft (collinear-usoft) effects is the purpose of 
this p^er. For recent a ppli cations of SCET in hard scattering processes and B-decays see 

Refs. BIlBHQEililEi- 

Since our focus is on mixed collinear-usoft interactions, we consider collinear quark helds 
collinear gluon helds usoft heavy quark helds usoft light quark helds and 
usoft gluons (We follow the notation in Refs. but for simplicity will often suppress 

the momentum label p on the collinear helds.) These degrees of freedom can interact in a 
local manner in Lagrangians and currents. This is in contrast with collinear-soft couplings 
whose interactions are mediated by ohshell huctuations 0], and appear in external operators. 
We comment on collinear-soft interactions at the end of the paper. 

The derivation of the leading order collinear quark and gluon Lagrangians and 
can be found in Ref. j3,|^, and a description of the gauge symmetries of SCET can be found 
in Refs. j3,0]. For details on power counting we refer to Ref. 13 • The heavy-to-light currents 
at large energy, J^i, were derived to leading order in Ref. p, including one-loop matching 
for all the Wilson coefficients. The running of these Wilson coefficients was considered in 
Refs. IJQ. 

In the context of the SCET, power suppressed corrections were hrst considered in Ref. 113. 
and the 0{X) suppressed currents Jm and collinear quark Lagrangians were derived. The 
authors showed that a reparameterization invariance (RPI) uniquely hxes the Wilson coeffi¬ 
cients of their sub leading currents and Lagrangian in terms of the leading order coefficients.^ 
In Ref. [ll| the RPI of SCET was extended to the most general three classes (I,II,III), and 
the multipole expansion of the collinear quark Lagrangian was treated to higher orders in A 
and were shown not to receive anomalous dimensions. In Ref. ^ the presence of additional 
0{X) heavy-to-light currents was pointed out that were missing in Ref. 3 |. 


The study of power corrections in SCET was continued in Ref. 3| and several important 
results were obtained for mixed usoft-collinear operators. In particular the mixed usoft- 
collinear quark Lagrangian C^g was hrst considered and was derived to (P(A, A^) working at 
tree level, but to all orders in attachments of ~ A° gluon helds. In a similar fashion 

heavy-to-light currents were derived to (P(A^), and linear combinations of currents that are 


^ A similar application of Lorentz invariance was used to derive constraints on the form of higher-twist 
contributions to structure functions in deep inelastic scattering in [l^ . For this case, invariance under 
changes in the light-cone vector fig was used to derive constraints on matrix elements (p|T'(/;(0)r'b(Ah^)|p). 


2 













invariant under the three types of RPI were identihed. It was also shown that the operators 
in C^q are not renormalized based on an analysis of arbitrary A^-loop diagrams in the hard 
region of QCD. The mixed usoft-collinear quark Lagrangian C^q was extended to a gauge 
invariant form with covariant derivatives in Ref. 


14 


The purpose of the present paper is to answer some open questions regarding our knowl¬ 
edge of the power suppressed usoft-collinear Lagrangian and heavy-to-light currents. This 
includes the number of Jhi currents at 0{X), since even at tree level the full reduction of 
Dirac structures has not yet been implemented. For both Jhi and C^q we also construct a 
complete basis which is valid for matching at any order in as, and therefore includes all 
operators that can be induced by radiative corrections or operator mixing. We work in the 
most general possible frame throughout (eg. allowing v± ^ 0, v-n ^ 1), and consider all the 
restrictions from RPI including the transformation of Wilson coefficients. Finally, we include 
the mixed usoft-collinear pure glue Lagrangian beyond LO (which follows from an extension 
of work in Refs, ji, ^). The above results are obtained by considering the full implications 
of RPI, and including all possible operators allowed from collinear gauge invariance, power 
counting, and the reduction of Dirac structures from the effective theory helds.^ 

For the heavy-to-light currents at 0{X) an important result we hnd is a new type of 
“three-body” currents, which have not been previously considered in the literature.^ In 
Refs. 0 , 01 the attention was restricted to SCET operators of two-body type J = 

... W){h„), where the two products in parenthesis are collinear gauge invariant, and the 
ellipses denote combinations of collinear derivatives. Beyond tree level but at the same order 
in A, we hnd that three-body structures can appear for some of the currents, having the form 
J = . hF)(lFl... W){h^) with three collinear gauge invariant factors. We show the RPI 

can be used to determine for which currents this happens. We also show that RPI greatly 
restricts the form of the three-body operators, so that they always involve a collinear gluon 
held strength. The two-body operators have hard Wilson coefficients which are functions of a 
single parameter while the new three-body operators have two parameter coefficients 

C(a;i,a; 2 ). Analogous three-body structures could appear in the usoft-collinear Lagrangian 
C^q at higher orders in perturbation theory, however using constraints from symmetries of 
SCET we prove that this does not occur. 

Our results are relevant to the study of decay channels for B mesons which involve 
energetic hadrons in the hnal state. For instance, the results derived in this paper are 
necessary ingredients in the factorization formula for heavy-to-light form factors proven in 
Ref. Q (for earlier work on factorization in heavy-to-light form factors see Ref. 00, 


^ Note that in deriving the complete basis for Jhi we restrict ourselves to 0(X) which is one order less than 
the order to which the tree-level matching results are known from Ref. [ijj. We treat £^q to O(X^), and 
give a detailed account of how t he gaug e invariant form in Ref. |lf| was derived. In cases where our results 
are restricted to those in Refs. |l0|,ll3 we find agreement, as discussed in more detail in the body of the 
paper. The results derived here are sufficient for the proof of a factorization theorem for heavy-to-light 
form factors to all orders in ag and leading order in 1 /Q E3 ■ 

^ In the final stages of this paper, Ref. [13 appeared where soft-collinear light-to-light currents are consid¬ 
ered. Although different from the usoft-collinear heavy-to-light case studied here, we note that 3-body 
currents were also found. Further remarks are left to a note added at the end. 
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and for results from QCD sum rules see Refs. Q) . The factorization theorem is valid to all 
orders in ag and leading order in 1/Q, Q = {rnsjE}, and separates contributions from the 
scales rsj Q\ ~ QA, and ~ A^, where A is a hadronic scale. It states that a gen eric 
form factor can be split into two types of contributions F = f^{Q) + f^^{Q) where |l5| 


f^{Q) 

r^(Q) 


pi pi poo 

Nq dz dx dr^T{z,Q,po) 

Jo Jo Jo 

X J{z, X, r+, Q, /io, p)(1>m{x, p)0B(r+, p), 


( 1 ) 


( 2 ) 


No = /b/m and the two terms both scale as This scaling is model 

independent and is in agreement with that derived from QCD Sum Rules [ 2 ^ . In Eq.__ 
and (j)B = (t>% are standard nonperturbative light-cone distribution amplitudes, c.f. 



The hard coefficients Ck and T can be calculated in an expansion in as{Q) and are simply 
related to the Wilson coefficients of the C>(A°, A^) current operators Jhi- The jet function 
J is dominated by momenta QA. If we wish to expand in as{y/QA) then using the 
techniques developed in Ref. J is calculable in terms of time-ordered products of the 
SCETi operators Jhi and C^g that we study here. At tree-level [ie. 0{aa{\/QAY agiQY)] 
one finds that J contains a 6{z — x), and in ratios of form factors the results for then 
agree with terms computed in Ref. |l2|. The z dependence first shows up at 0{ag) as does 
possible dependence on 0^. However, as we show in section ED it is possible to absorb the 
terms into a redefinition of the to all orders in perturbation theory. 

The factorization formula provides a clean separation of the “soft” non-factorizable (NF) 
contributions and “hard” factorizable (F) terms without double counting. It also gives us a 
procedure to systematically improve the predictions to any order in perturbation theory at 
leading order in 1/Q. The value of T and Ck can depend on which heavy-to-light process 
we consider, whereas 0 m and 0| are universal functions. The (kS are also universal since 
only a C^{E) appears for decays to pseudoscalars M, and a Cx{E) and C\i^{E) appear if M 
is a vector meson. The jet functions J are common among certain classes of form factors 
and also do not depend on the precise state (eg. n or rj). The terms satisfy the so- 


called large energy form factor relations 
these terms as “soft” contributions \A llS , 


as expected from the prior loose definitions of 
|23[.^ Note that we have not bothered to separate 
~ QA and p^ ~ A^ fluctuations in the functions, since it is not clearly beneficial 
phenomenologically. The factorization theorem does tell us that ~ (A/Q)^/^, however it 
does not distinguish between factors of irib and E in this Q~^l’^. It also does not numerically 
favor the or term, for instance it is possible that the leading as{y/QA) in J is 
compensated by an analogous factor in (//f. 

We start in Section 11 by reviewing the general constraints imposed on SCET operators 
following from collinear gauge invariance, spin structure reduction, and reparameterization 
invariance. In Section 111 we study the implications of these predictions for the subleading 


These relations were first derived in Ref. [ 2 ^ using LEET 0. However for studying energetic hadrons 
with QCD the LEET framework is known to be inconsistent P, Eil. for instance it does not bind an 
energetic quark-antiquark pair into a meson in heavy-to-light decays |2fi| . 
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usoft-collinear Lagrangian Cue- In Section IV we present detailed results for SCET currents. 
Using the example of the scalar current as the pedagogical example, we demonstrate the 
construction of the complete basis of 0(A) operators contributing to the weak currents, 
which closes under RPI transformations. Explicit results are then also derived for the 
pseudo-scalar, vector, axial-vector, and tensor heavy-to-light currents to 0(A). In Section 
V we summarize the one-loop matching results for the currents, give explicit results for C^g 
Feynman rules, and discuss the basis of currents in the particular frame v_[_ = 0, n-v = 1. 


II. OPERATOR CONSTRAINTS IN SCET 


In this section we briefly review the symmetries and structure of SCET which will be 
important for our construction of operators. We refer to Refs. Ed| for more 

details. 

SCET includes infrared degrees of freedom corresponding to the relevant low energy scales 
in the problem. These are typically those with momentum that are collinear Q(A2,1,A), 
soft p^ ~ Q(A, A, A), or ultrasoft (usoft) ~ Q(A^, A^, A^), where the components here are 
in a light-cone basis (-1-,—,T). Each type of mode has effective theory quark and gluon 
helds, which are then organized into operators with a well-dehned power counting in A. It 
is convenient to introduce light-cone unit vectors satisfying ri^ = = 0, n-n = 2, 

in terms of which a vector has components p^ = {n ■ p,n ■ p,p^)- The couplings of the helds 
are described by an effective Lagrangian, while the couplings to external sources appear 
as additional operators or currents. Both the Lagrangian and currents are constructed 
such that they include constraints from power counting, spin symmetries, and collinear and 
(u)soft gauge invariance. 

The soft-collinear effective theory also contains a kinematical reparameterization invari¬ 
ance symmetry. Lorentz invariance is broken by introducing the vectors n and ft, but 
is restored order by order in A, by requiring invariance of operators under a simultane¬ 
ous change in n and/or n and compensating changes in the effective theory helds. This 
reparameterization invariance (RPI) symmetry of SCET was hrst considered in Ref. H. 
and was then extended to the most general three classes (I,II,III) of allowed transforma¬ 
tions in Ref. 
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The three types are dehned by the inhnitesimal change they induce 
on the light-cone unit vectors: type-I {n ^ n + A_i_), type-II (fi —>• n -f e_i_), and type-III 
(u —> (1 -|- a)n,n —>■ (1 — a)n). Here a ~ ~ A°, while ~ A. It is the analog of the 

reparameterization invariance of heavy quark ehective theory (HQET) under changes in the 
heavy quark velocity v ji^, where = 1. We will use HQET for heavy quark helds Q 
The restrictions we consider for hnding the most general set of power suppressed gauge 
invariant operators are: 


® The nature of Lorentz symmetries on the light-cone are well known . The new point in SCET |l(llll| 
is that for any collinear process these symmetries are realized in a way that leads to non-trivial restrictions 
both on operators at a given order in the power counting and between operators at different orders in the 
power counting. 
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i) Power counting and gauge invariance which determine what basic building blocks are 
allowed at the order we are considering. 

ii) What auxiliary vectors are available (such as n, ft, n, ...), which can be used to 
construct the most general set of allowed scalars/tensors/Dirac structures. 

hi) Eliminate operators which are redundant by integration by parts, or equations of 
motion. 

iv) Impose type-III reparameterization invariance. If a non-trivial invariant can be formed 
with the label operators, such as (n-vV), then include Wilson coefficients that depends 
on these quantities. 

v) Impose all constraints from type-I and type-II reparameterization invariance. 


To impose the hve constraints we start by writing minimal sets of independent operators 
compatible with the general principles in i), ii), and iii). We then require RPI invariance 
order by order in the A power counting. To do this we found it useful to split the RPI 

(A°) 

transformations into two categories, those that act within the order we are considering 5) , 

with j =1,II,III, and those which connect operators to one higher order At leading 
order the type-II and type-III o) transformations already provide non-trivial constraints 


on the allowed form of operators. In contrast the transformations allow us to derive 

erators at different 
iOl , where we note 


relations valid to all orders in as between the Wilson coefficients of or 
orders in A, These relations are similar to the case of RPI in HQET j28|, 
in particular the important relations derived for coefficients of subleading heavy-to-heavy 
currents in Ref. [s^. We start by summarizing restrictions that follow from collinear gauge 
invariance and power counting in section III Al spin structure reductions in section mu and 
RPI in section Em 


To separate the momentum scales we follow Ref. and use collinear quark helds ^n,p{x) 
(and gluon helds A!^p{x)) which have momentum labels p for the large components of the 
collinear momenta, and residual coordinates l/A^I . Thus, all derivatives on collinear 
helds are the same size as derivatives on usoft helds, d'^ ~ A^. This setup implements 
the multipole expansion in momentum space. Note that our analysis of power corrections 


dihers from Ref. 113 in two ways, the hrst being that in Ref. 113 the momentum scales 
were separated by performing the multipole expansion in position space, which however 
leads to an equivalent formulation. We do hnd that concise results for the power suppressed 
corrections are obtained with the momentum space version. Secondly, we derive our basis 
of operators and implement all symmetry constraints working order by order in the power 
counting, rather than constructing invariants and then expanding in A. This made it simpler 
to derive a complete gauge invariant basis at the desired order while working in a general 
frame. 
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A. Power Counting and Gauge Invariance 


The SCET is derived from QCD by integrating ont flnctuations with ^ where 

in typical processes A = (Aqcd/Q)^ with fc = lorfc = l/2. Infrared fluctuations are then 
described by effective theory fields. A gauge invariant power counting for fields can be fixed 
by demanding that the kinetic terms in the action are order A°. For the collinear fields 
this gives ~ A*^ for the quarks, and {n■An,n■An,A^^^) (A^,A° A) for the collinear 

gluons, hv ^ q ^ y? for usoft quarks, and ~ A^ for usoft gluons |l|, |^. Derivatives on 
these helds count as ~ A^. The larger collinear momenta are picked out by introducing 
label operators "P ~ A° and ~ A j3|. For example P in,p = {n-p)^n,p- For notational 
convenience we define collinear covariant derivatives 


in-Dc = V + gn-An, = V^ +gA^j^^ , (3) 

and ultrasoft covariant derivatives 

in-Dus=in-d + gn-Aus, . (4) 

For the components, it is only the combination 

in-D = in-d + gn-An + gn-Aus , (5) 


that ever appears. In general a derivative without a subscript involves the sum of the 
collinear and usoft pieces, and it is this combination which is RPI invari¬ 

ant HD (implying that the anomalous dimensions of terms that appear in the multipole 
expansion are related). 

Integrating out the offshell fluctuations builds up a collinear Wilson line, W, built out of 
collinear gluon fields which are not suppressed in the power counting 


W = [ 5^ exp £ n ■ A^Jx) ) 

perms 


( 6 ) 


where the label operators only act on helds inside the square brackets. Up to the important 
fact that W has been multipole expanded, it is the Fourier transform of a standard position 
space Wilson line, W(—oo,x). Factors of lU ~ A° can be included in operators without 
changing the order in the power counting. However, their location is restricted by collinear 
gauge transformations, Uc, under which W —>■ UcW j^. Since P ~ A° in the power counting 
the hard Wilson coefficients can be arbitrary functions of the momentum or momenta, a;*, 
picked out by these operator, Ouji^p) [^. These coefficients can be computed by matching 
with QCD at the hard scale /i ~ Q and running with the renormalization group. 

If we consider a general Wilson coefficient and operator C then the covariant deriva¬ 
tive 


in - Dc = WPW\ (7) 

so it is always possible to put all the Wilson lines in O and the dependence on the momenta 
picked out by P into C. We will find it convenient to use the notation 

{^nWU = [UW6{uj,-n-vP^)] , 

6{u;2-n-vP^)] , ( 8 ) 
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where again the label operators do not act ontside the sqnare brackets. The factor of 
n-v is inclnded next to V to make it a type-III RPI invariant. Thus, the momentum 
labels oji do not transform under RPI. The products of helds in Eq. (jS)) are color singlets 
under the collinear gauge symmetry, so the momentum labels uji are gauge invariant. These 
products still transform under an usoft gauge transformation Uus as (^nkP) —>■ (■CnfP) 
and {WW^^W) U^s{WW^^^W)Ul,. We will elaborate on how RPI affects Wilson 

coefficients in SCET in subsection III (Jl below. 

For Lagrangians and currents where the variable is not available we can not make use 
of the dehnitions in Eq. (jH)). It is still convenient to make use of a similar notation: 

(CnkP)., = [lnW5{z^-V^)\, 

= [W^D^^W 5{z2-V^)\, (9) 

where we use the variables Zi rather than a;*. Under a type-III transformation the Zi trans¬ 
form like n so the delta function is homogeneous (and compensated by an integration measure 
dzi). 

Using the scalings for helds and derivatives the power 
A^, can be determined entirely from its operators using 

S = 4+5^(«:-4)[Vf+ ■ (W) 

k 

Here count the number of order operators which have collinear helds, soft 

helds, both, or neither respectively. For any operator the power of k is derived by adding 
up the powers of A in its components, so for instance ~ A^, counts as = 1. 

Since the operators are gauge invariant so is their value of k and also the power counting 
of any diagram using the result for <5 in Eq. (uni). In this paper we focus on operators with 
ys = ysc ^ 0. 

We have also found it convenient to dehne additional pure gluon operators. In particular 
we will use the purely collinear held strength 

( 11 ) 

We will also make use of the mixed tensors 


counting for an arbitrary diagram, 

31 


1 n rr-t . A/f 


In fact the operators Ijlli and n ■ M, together with ig n-M = [in-D'^, can be combined 

into a single object closed under usoft Lorentz transformations, which transforms in the 
desired way under the collinear and usoft gauge symmetries 

igljld= + ^gn-An ]. (13) 


Finally the following results for manipulating covariant derivatives on Wilson lines also prove 
to be useful 


(wUpjw) 
(W^ ipjw) 


[wt ipiw] -f[ = 

[Wt ipjW] +f± = 


-1 

-V 
-1 

-V 


W^igp^W 




WHgp^W +fx, 


(14) 
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B. Reduction in Spin Strnctures 


Collinear quarks and heavy usoft quarks have spinors with only two non-zero components. 
In four component notation this is encoded in projection formulae for the helds, 


Priori Cn 




( 15 ) 


where and = {l + fi)/2. We also dehne the orthogonal projector P^ = 

)/4 where = 1. A quark bilinear with a heavy ultrasoft quark and light collinear 

quark therefore only has four possible non-trivial Dirac structures. On the other hand if the 
heavy ultrasoft quark is replaced by a massless ultrasoft quark which has a four component 
spinor then their are eight possible Dirac structures. When generating operators we should 
be careful not to include redundant Dirac structures. Therefore, it is convenient to have a 
canonical basis which we can project results onto to check their interdependence. For this 
purpose we choose the basis 


1 Fi 

^ 71^2 Qus ) ^2 


2’~’ 

1 * 7= 

’2’^ 



7±, 



(16) 


Any general Dirac structure can be projected onto a linear combination of terms in these 
basis with the help of the following formulae 


and 


(^^^F hqj hy , 

2 L2 


ri = -tr ^PnFP^ 


t^tT[^^,PnTPy 


+ tr 


'y-^P-TP 

I u V 


(Jus ^ 71^2 (Jus •) 

U = itij^r 

L 8 

I A* f 

+77 tr 


^ \'4 


H— tr 
2 


^F 
4 . 

^7( 


+ 75 tr 

tr 


#75 #5 , (#5 


( 17 ) 


( 18 ) 


The number of independent structures is quite logical, for ^rJ'hy each held is determined 
by two-component spinors and there are 2x2 = 4 terms in the basis. For ^rJ^Qus only the 
collinear spinor has two-components and there are 2x4 = 8 terms in our basis. Our choice 
of basis in Eq. (ED differs from Ref. where the choice Fi = { 1 , 75 , 77 } was used, and 
calculations were given in a frame where v-n = 1 and v± = 0. When is kept arbitrary we 
have found the basis in Eq. (Hi is more convenient since it retains its orthonormality in an 
arbitrary frame. 

The projections formulae in Eq. ED can be used to reduce the possible Dirac structures 
in constructing a complete basis of operators. It is convenient to dehne r7 = 
and u7 = — n-vn^/2 — n-vn^/2, since then ^7; ''"±1 form a complete vector basis. 
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To reduce the Dirac structures we can use relations such as 


^_L = 1 


le 


fiuafS 


n-v 

I 

ilh 




^_L = — 2n-u , 


7± = 7m 


n-v 




r^75 = - 2 u^ + 27 ^ + ^7^77,-u , 


^’"±75 = — 4n-u7^ — 2 ^ 7 ^ , 


naV/3 - {v^'yu-Vy^^) - 


n-v 




( 19 ) 


where the = indicates that these are only true between the fields in Eq. jni). (The complete 
set of relations is rather lengthy and is not shown.) The relations in Eq. (unD allow the 
structures on the left to be traded for those on the right (with more than one iteration 
in some cases). Using the projection formulae it is straightforward to show that the most 
general Dirac structure possible for the LO scalar currents are {1,'^}, while the vector and 
axial-vector currents have the basis shown in Eq. (HI, and the tensor currents depend on 
the basis in Eq. (El 


C. Reparameterization Invariance 


The decomposition into collinear fields requires introducing two light-like vectors n and fi, 
such that = 0 and n-n = 2. These vectors break five of the six Lorentz generators. 

This part of the Lorentz symmetry is restored order by order in the power counting by 
requiring invariance under reparameterization transformations on n and n 
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(I) 



( 11 ) 


nf, 

^ 


(III) 


n^ ^ {1 + a) n^ 

^ (1 - a) n^ 


( 20 ) 


where A-*- ~ A, while ~ a ~ A°. In general one has two options for constructing 
RPI invariants: i) construct operators out of completely RPI invariant quantities and then 
expand these in powers of A, ii) construct operators order by order in A and transform them 
to see what linear combinations are invariant, and which operators are ruled out. In this 
paper we will adopt approach ii), since starting with the most general gauge invariant sets 
and then reducing them allows us to be confident that we do not miss operators that could 
arise at any order in perturbation theory. 

For our purposes it is convenient to divide the RPI transformations into two subsets, 

(A®) (A*^) 

those which include terms within the same order in A denoted by , and hfjj , and 

those which cause order A suppressed transformations denoted by and All type- 
III transformations act within the same order in A and it is easy to construct invariants 
under type-III. We simply need to have the same number of n’s (h’s) in the numerator and 
denominator, or have products of n times n. The transformations of type-I and type-II are 
more involved. From Ref. the transformations that have terms of the same order in A 
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are 




hi 


n-D ^ n-D + Dj 




A'; 

—— n-D , 
2 


( 21 ) 


n- 




n'^ 


n^ + £((, - 7± ^ 7± - , 


and the transformations that start one power down in A inclnde 






+ A); , 




u u 

'll ^ hi 


A') 


n 


i - —, 


r{-^) 

hi 


2 ' 2 
ft-Zl —>■ h-Zl + £^-iA_L , i7(( 

1 




■Cn 'Cn ( 1 + Pl- 


n 


■D 2 


W 


n-D 


e^-D^]W 


( 22 ) 


( 23 ) 


Here W is the RPI completed W, and is the Fonrier transform with respect to y of a position 
space Wilson line involving [fi-An + n- Aus){sn + x) taken from s = —oo to y Q . When 
expanded in A, W = IF + (P(A^), where IF involves only the n-An held as in Eq. (jHl). 

If we start by considering LO operators then they must be invariant under the 
transformations in Eq. m all by themselves. The transformations of the LO terms 
connect them to NLO operator’s transformations. Since in the collinear sector only 
terms exist this pattern repeats at all higher orders in the power counting. Note that 
here we will not need to consider HQET RPI under the velocity v^. Since the transformation 
^ + A(( where A(( rsj Aqcd/Q rsj A^, this type of RPI only needs to be taken into 

account at one-higher order than the order we are working. The combined SCET and HQET 
RPI transformations were used in the (P(A^) analysis of Jhi in Ref. Q. 

Finally we consider a new feature of RPI in SCET, namely how Wilson coefficients are 
affected by reparameterization invariance. Our analysis is similar in spirit to Ref. js^, 
where heavy-to-heavy HQET currents with coefficients depending on the change in velocity, 
C{v-v'), were analyzed. If we adopt the view of building invariants at all orders in A then 
the coefficients in SCET must also be functions of invariants, such as operators like 


^nC{-iD-V)Tn^ 


( 24 ) 


where Hv are invariants including the quark helds hy, and is the RPI version of 
the velocity ji^. When expanded in A the leading term involving the covariant derivative 
in C can be traded for W and V using Eq. 0,c( -in-vn-Dc) = WC{n-vV^)W'^. Here we 
will use the opposite but equivalent arrangement of starting with a current that is leading 
order in A, 

^nWC{n-vV^)K, ( 25 ) 


and then determining how both the operators and coefficient transform under RPI. We then 
determine which structures are required at one higher order in A to cancel this change, and 
which allowed higher order structures are left unconstrained. 
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III. COLLINEAR-ULTRASOFT LAGRANGIAN 


In this section we discuss the mixed ultrasoft-collinear Lagrangians to These 

actions are power suppressed [^, and start at 0{X) 113 . In section mm we consider the 
derivation from integrating out components of the full theory field, which gives a tree level 
derivation of the action (for further explanation of this approach see Refs. BBE3 ). In 
Ref. 0 this procedure was used to derive a form for the mixed ultrasoft-collinear quark 
Lagrangian, but a manifestly gauge invariant form was not determined. In Ref. the 
analysis was extended to give manifestly gauge invariant operators in terms of covariant 
derivatives. In section nnn we review the details of how a derivation of a gauge invariant 
form of the action was carried out in Ref. 0 where the result is purely in terms of field 
strengths. 

However, since the analysis in section IIII Al is only valid at tree level, it misses i) non¬ 
trivial Wilson coefficients in the tree level operators, and ii) new operators whose coefficients 
can have zero tree-level matching. In Ref. point i) was addressed and it was shown 
diagrammatically that no non-trivial Wilson coefficients are generated. However, point ii) 
has not yet been addressed, so additional operators could still be induced by matching at 
some higher order in perturbation theory. In section ITHRI we show that both points i) and ii) 
can be simultaneously solved by using the full set of symmetries of SCET when constructing 
operators. We also extend the derivation to the mixed usoft-collinear pure gluon sector. 


A. Matching for Cue at tree level, but all orders in n-An gluons 


In this subsection we discuss in detail the matching calculation for the mixed usoft- 
collinear quark Lagrangian 1^. The part of our discussion from Eq. (OH]) to PI|) follows 
Ref. US , but with our momentum space notation. We start with the action C = 'ipip'ip and 
decompose it with SCET fields 


- fL 

C' T T ^nQ-^-c Qus T Qus QP-c T Qus QP-c T (Jus 'i'P^us (Ju 


+ 


_ 4 

T ^n9-P-c(Ju 


( 26 ) 


where the D is usoft plus collinear, D^. is purely collinear, and collinear momentum conserva¬ 
tion has been enforced. Varying with respect to gives an equation of motion to eliminate 
this field from the term in square brackets 




1 y/i 


in-D 2 . 


T 9-P-n (Ju 


^n = 


(JusdP-n 1. 


^ 1 


in- D 


( 27 ) 


Plugging this into Eq. (EHl) and expanding we find that the two collinear quark terms exactly 


reproduce terms in the gauge invariant multipole expanded action in Ref. 
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Using 


® Note that in QED the C(A) pure collinear quark Lagrangian can be written in terms of With 

the momentum space multipole expansion this manipulation is not necessary to achieve a gauge invariant 
result 
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Eq. (jlOj) the terms with two ultrasoft quarks and > 2 collinear gluons hrst show up at S 
ie. C>(A3) , and are therefore neglected. The mixed usoft-collinear quark terms are 




^nQ-^nQus j-^d^nQu 

2 m-U 


+ 


Q.US 9-^ri T Qus 94 


n . — 

in-D 




Taking Eq. (EHI) and expanding to second order in A gives 


= 3 , 

(28) 


= in[94\ - 9n-A^qus + h.c. , (29) 

^fq = iX ( 9 n-A^ + 941)(1 us - 9 n-A'' q^s + h.c., 

where the superscripts denote the power suppression in A and A^. In each of the last terms in 
1 — IT = l/(ih-iAc) qn-Ac, and in the hrst term of we write g4\ = 
and then q^s = [IT + (1 — IT)]g„s- Thus, 

= ln{^P^_W -fi)qus + h.c., 

4? = en|(^/n-Al=+*.^i^^(*.^I-n))[IT+(l-IT)]g„, + en*|)r(W^-l)?«. + h.c. 

= ipl)wq^s - n Qus 

+ zK.*)(IT-l) + h.c.. (30) 

In manipulating we used the fact that integration by parts is allowed on the (1 — W)qus 
term and we can then use the equation of motion for the collinear quark to give a term 
{-in-D us)4 ~ IT) which we collected with the — 1) term. The result in Eq. (pn|l 

agrees with Ref. [l3| , up to the fact that we performed the multipole expansion in momentum 
space. 

In Eq. (j3()|l we did not drop the ^±9 = ^ terms since we want to make explicit the fact that 
it is the combination {ip]_W — = [i.^IT] which starts with at least one-collinear gluon. 

Written this way it appears that our is not collinear gauge invariant. In the transformed 
result the non-invariant term cancels if we use 'f’i_qus = 0, but then it is not explicit that 
the operator starts with one-collinear gluon, so has either one or the other explicit. 
For Eq. (HirH) still involves the gluon field A^ so the gauge invariance of this expression 
is not at all clear. However, the above considerations indicates that it should be possible 
to write all the terms in Eq. dsni) in terms of gluon held strengths, and thereby achieve a 
manifestly gauge invariant action that starts with one-collinear gluon. This derivation was 
carried out in Ref. jl^, but no details of the calculation were described there. These details 
are described below in Eqs. m through dsni). 

To proceed we note that using Eq. (ITTll . igB^j_W = [in-Dc,iDd]W = in-DciD^W — 
iD^WV. Making similar manipulations for n-M and Ifii we can write 


igp^W = in-D, {il/pw - fP) - {Hjpw - 'fA}P + {gn ■ Apf^ . 
ign-MW = in-Df-in-DW — Pin-D^s — {in-DW — in-D^sjP , 
igpI^AV = in-D^ippW - Pipp + {ipp - ippW}P. (31) 
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Now we take purely usoft fields on the right, and divide on the left by in-Dc- In Eq. m 
the terms in curly brackets start at one-collinear gluon, so even in the presence of 1/{in-Dc) 
these terms are non-singular and can safely be dropped using that the label operators give 
zero on the usoft held. This gives 

= {i]ptW-'t)q^,, 

= {in-DW - in-Dus)qus - {W - l)in-Dusqus , 

= - (IE - 1) ip^lq^s . (32) 

These expressions allow us to write covariant derivatives acting on Wilson lines in terms of 
held strengths. 

Using Eq. for in Eq. (IHn|l . we arrive at the hnal result 

4? = ig^Wqus + h.c.. (33) 

TjTIj * -L'c 

This form is particularly nice since it is explicitly collinear and usoft gauge invariant and 
furthermore explicitly starts at one-collinear gluon due to the To see the gauge invari¬ 
ance note that under a collinear gauge transformation Uc we have Udm W —> f 41 U, 

^ Uc$lUI and {n-Dc)-^ ^ Uc (h-D,)"^ U] so all factors of Uc cancel. Under an 
ultrasoft gauge transformation Uu we have —> Uu^n, hU — U^WUl, ]/^±_ Uup±Ul, 

{n-Dc)~^ Uu {n-Dc)~^ Ul, and qus Uuqus so all factors of Uu also cancel. In Fig. El in 
section rV (II we show the one and two gluon Feynman rules that follow from Cue in Eq. O- 
A non-trivial check on our manipulations is that the same Feynman rules can be obtained 
from Eq. (HB by using the free equations of motion. 

We now proceed to further simplify C^^^ in Eq. Using Eq. leaves 

4'^ = qus + inl^]plJT^^^g$lW qus 

in-Dc 2 [in-Dc) 

+ln{W - l)ipus qus + h.c., (34) 

where in the hrst and last terms we used the fact that = 0 to write a full 7 ^ in and 
Ipus- For the last term in Eq. we can now use the equation of motion for the usoft quark 
held to give our hnal result 

4? = igm qus + WplW qus + h.c. (35) 

in-Dc 2 [in-Dc) 

Again in this form the action is explicitly collinear and usoft gauge invariant and fur¬ 
thermore explicitly starts at one-collinear gluon due to the held strength i?’s and M’s. In 
the way we have written the result it is invariant under usoft Lorentz transformations on 
which separately rotate 7 ^ and 

Finally we note that the mixed usoft-collinear quark actions in Eqs. (El and dSl proved 

to be important for the proof of a factorization formula for heavy-to-light decays in Ref. jl^ . 

fl 2) 

In the next section we analyze the most general possible basis for ’ beyond tree level, 
which follow purely from symmetry considerations and also discuss power suppressed terms 
in the collinear gluon action. 


in-Dc 

._4 ign-MWqus 
in-Dc 

.-4 igpIiWqus 
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B. Most General Basis for C 


UC 


The ultrasoft-collinear quark Lagrangian can be expanded in a power series in the param¬ 
eter A. It is not possible to construct an invariant operator that is dimension-4 and order 
A°. Therefore, we have the series 

+ • • • • (36) 

Since this is a Lagrangian, insertions of these operators do not inject momentum, and we 
are free to integrate by parts as long as we are careful not to generate singular terms. 

To construct the most general quark action we can use a single collinear quark held 

~ A, an ultrasoft quark held Qus A^, and a A. These factors give a dimension 

4 operator, and from the power counting formula in Eq. m they give 5 = 1 which is the 
correct order for ■ To satisfy collinear gauge invariance without changing the order in the 
power counting we make use of the Wilson line W to write (^nhh) and Since 

the Lagrangian is a scalar we must dot the index ji into another vector. The possible Dirac 
structures are restricted by the fact that ir^^n = 0. They are also restricted by type-III RPI, 
for instance ^ 7 ^ is not invariant and therefore is ruled out (in the case of the heavy-to-light 
currents we can make use of the product n-v, so n-viji is allowed). Taking these constraints 
into account leaves 7 ^ as the only possibility. Thus, we have reduced Cue to the form 

iLW) p{V, p) {W\1P^W) qus + h.c., (37) 

where the coefficient p is dimensionless and ilp^ acts to the left or right. However, by type- 
III RPI invariance the coefficient p('P, p) can not be a function of V, leaving only p{p). Now 
p is a dimensionless function of the dimension-full parameter p and can only be equal to 
a constant (assuming no new dynamical scales like Ageo are generated by renormalizing 
c\i}'). Now lp±q = 0 since q carries no perpendicular momenta of order A, so using Eq. m 
we see that can be traded for a Ijij: operator. For we obtain the same 

Iji^ operator, plus = 0 ; which follows from the fact that and by momentum 

conservation {^nW), carries zero collinear T momentum. Fixing the constant p = 1 by tree 
level matching then leaves 

cf^ = ^nW^WUg^Wqus+ ig^Wqus + h.c. . (38) 

V tn-Dc 

In this form it is clear that the operator is collinear and usoft gauge invariant and generates 
terms with > 1 collinear gluon as required by momentum conservation. Finally, it is easy 
to show that its order A*^ type-I and type-II variations vanish, since from Eq. m only the 
transformation of must be considered and Sf' oc [in-Dc, in-Dd = 0 , while ^ 

which gives zero since = 0 . 

The above line of reasoning can be repeated at (P(A^). By power counting and gauge 
invariance we can now have one {WHn-DW) or two {W'^Dd^W) factors with derivatives 
to the left or right. Again we can use Eq. (jl4|l to simplify the covariant derivative terms. 
We must have an operator starting with one-collinear gluon, and again the iW'^DW) factor 
next to qus must be in square brackets and can be turned into a gluon held strength. Also 
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by type-III RPI the Wilson coefficients must again be numbers, except for operators with 
three or more invariant collinear products in where they can be functions of the ratio Zi 
of minus-momenta. Taking into account these constraints leaves three possible operators ^ 

= Pi {LW)(w^^^igI^w)qus + h.c. , (39) 

V m-Dc / 

cf = jdzidz, 

Cf = jdzidz,p3{^) (lnWU^(WHD^'‘WU~{w'^igB^w) q,. + h.c. , 


where (.. = [... 5{zi — P^)], and we’ve used our freedom to integrate by parts to make 

the perp covariant derivatives act to the right. Note that the overall n-p momentum is zero, 
so no Zi label is used on the l/ij: bracketed term (following the convention in Ref. [^). Again 
the presence or absence of factors of ifi are completely fixed by type-III RPI. Now consider 
the type-I and type-II RPI transformations. Computing the order A variations of Cu} and 
simplifying the resulting expressions gives 


c 


ilp. 


2 in-Dr 2 


m-D, 




(w-T>c)^ 




in-Dc 2 


[6f\gm-)Wqus + h.c. , (40) 

ign-M+^^^[6[pigIjliL)\Wqus + h.c. 

ZTZ * j—'c ^ 


The igljfiL terms appear since in the RPI transformations in Eq. it is the full 
which transforms. For the order variations of we hnd 


r(A*^) /»(2) ^ 

61 ’C \' = Pi 

r(A*^) /^(2) _ ^ 

^1 ~ Pl ^r, 


1 1a^-b^ + 


in-Dc2 

1 

in-Dr 2 


-^{Sp^igJ^liLpWqus + h.c., 

11L' Jr.y Q ^ 

ign-M + ^J—[5[f\g]^Pl\Wqus + h.c., 

%T\i * -L-^n ■' 


(41) 


^dp = - dzidz2 P2 ( 


2^1 

(- 

Vzi 


I 

4 V 

dnW).Aw^ip)iw) 


ip’rf = -ldzidzzP2(^) (LW)2Az2)^6{w'^^ujftw)q^ + h.c., 


'2 p\ in-D^ 


+ [dz,dz 2 P 2 (-) dnW),q{wdpiW),,i(w^^^igerB^w)qq,, + h.c., 

J \zi/ V\ in-Dc / 


\iD) ^(2) 


fZ2 


6P’C);> = -JdZidZ2pp^ 


{inW),,6{z2) 


r('^°)p(2) _ 


4 "P V in-D^ 
1 


igA -B^W]qus + h.c., 
1 


dz2dz2P2D) 


^ This assumes we have eliminated a possible four quark operator using the collinear gluon equations 
of motion g‘^{^nWT^^Wpn)l/V‘^{^nWT^i^qus) = ■ Dc) ign - MWqus + 2^nT^/{2{in- 

Dc)^} [iDr^,igBr^]Wqus- This conclusion is not changed if we consider the most general possible four- 
quark operators allowed by all the symmetry constraints. 
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Comparing Eqs. iOD and (ED we see that it is not possible to form an invariant involving 
£ 3 ^^ so p 3 = 0, while an invariant can be formed from Cue + Cf^ + C^ 2 ^ by taking pi = 1 
and P 2 {z 2 /zi) = 1. Since p 2 is independent of zijZ 2 the integrals over 2 : 1^2 can be performed. 
Therefore, we can write our hnal result for the hrst two orders in the usoft-collinear quark 
Lagrangian as 

4 ? = + h.c. , 

tn-Uc 

ig^ w qus + h.c., 

c = ^ 

These terms agree exactly with the result from tree level matching in Eqs. (El and (El- 
The analysis here shows that no other terms are induced by matching at any order in as- 
Next we proceed to analyze power suppressed terms in the collinear gluon action. Starting 
with the LO collinear gluon action |^, and making it RPI invariant with iD'^ = iD^ + 
gives 

. (43) 

It is straightforward to see that no other gauge invariant pure glue dimension-4 operators 
are possible. We could build a more general gauge invariant operator out of a string of m 
terms with m — 4 factors of l/P to make up the mass dimensions. However, 

type-III RPI then demands m — 4 factors of n^^ in the numerators which using n-DyV = 0 
collapses the operator to the case m = 4. Finally, since W transforms under type-II RPI 
as in Eq. El, but iD^ does not we hnd that these operators must have Wilson coefficients 
C{zi) that are independent of the Zi parameters. In this case all factors of W cancel out 
and we are left with Eq. El (after performing the Zi integrals and hxing the coefficient at 
tree level). Expanding Eq. (HI we see that the order A and suppressed terms are 

^ I®", J } . (44) 

■CS' = ^ I®". ®i''] K. ®t J} + y I®™. ®t1 [®i. ®i]} 

where = D'^ + n^^n-Dus/‘2. 

In Ref. 0] the gauge hxing terms in the LO gluon action were given in a general covariant 
gauge. We do not bother to consider the possibility of other leading order gauge hxing 
terms since we have some residual freedom to choose these terms however we like. In an 
RPI invariant form the terms from Ref. yi are 


Ccg = 2tr 


Wf,, iDf^,Cn I^ tr I , H);] 


(45) 
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and the subleading terms in their expansion are 


/•(i) 

^cg 


^cg 


= 2 tr< c 




, Cn 




, Cn 


+ 2 tr| Cn 

iD:l;;^,Cn | + tr| c„ in-D'^% in-D,Cr, 
in • Dus 5 Cn ^ 

+ ^ tr { [Wl^^ , [Wll , A^’'] I + i tr I [in-D^^ ,n-An] [iD^ , A";^ 


2 

+ - 
a 

2 tr<i c. 


(46) 


+ tr| Cn 
1 


jr)us 

in-D, 


(47) 


IV. MOST GENERAL BASIS FOR HEAVY-TO-LIGHT CURRENTS 

In this section we give our derivation of the most general basis of heavy-to-light currents 
at 0{X). The scalar current is given in great detail, and forms the basis of the analysis for 
the other Dirac structures. Expanding the heavy-to-light currents in powers of A we write 

Jm = + ..., (48) 

for the LO currents (J^*^^), and NLO currents The superscript denotes whether the 

current is scalar (d = s), pseudo-scalar (d = p), vector (d = n), axial-vector (d = a), or 
tensor (d = t). For the preliminary basis where only constraints from gauge invariance, 
power counting, and type-III RPI invariance are imposed we use a calligraphic notation 
and then switch to Roman for the hnal basis that is invariant under all 
the type-I and type-II constraints. We will also make use of a convolution notation 

= [du (a;), (49) 

J \mh mbJ 

where contains helds and operators with the notation in Eq. (jH)) and the Wilson 

coefficients ^ p/m^) are numerical functions of the convolution parameter (where Cj = 

uj/mb). 


A. Scalar Currents 


From gauge invariance and power counting the most general leading order heavy-to-light 
current has the form j^. For n_L 7 ^ 0 the most general allowed scalar spin structures 

from section mu are then T = {1,'^}. Type-III RPI demands that the ^ is accompanied by 
either a.n-v or a 1/n-v. Thus after imposing constraints i)-iv) of section ITTl we are left with 
the possible leading order currents 


= 


(^) 

= Cl 




= 


(s) 
= C2 


AVn-v^h^ 


= cT CnhP 


(*) 




2n-v 


~hq) 


(50) 


where q are dimensionless Wilson coefficients. With type-III RPI invariance the q can only 
depend on the combination {n-v V), the b-quark mass mb, and the renormalization scale p. 
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Now consider the order type-I and II RPI transformations in Eq. m- Since none of the 

(A°) 

operators in Eq. ()5()|1 involve qnantities that have dj ^ transformations they are invariant 
nnder type-I at this order. However, nnder the analogons type-II transformations 


= 0 , 




^nlV 


i 


2n-v 2{n-vY 


e^-v]h^. 


( 61 ) 


Thns, it is not possible to form an invariant involving the cnrrents only the cnrrent 

is allowed. Therefore, we can rewrite onr hnal result for the most general leading order 
scalar current as 


t(^) _ r^(^) 

"'O ~ *-'0 


-n-vV /i 

mb ’ mb 




( 52 ) 


Since the Wilson coefficient is dimensionless it can only be a function of the dimensionless 
ratios of parameters as shown. The minus sign in the hrst variable is included so that V gives 
the total outgoing momentum of Switching to the convolution notation in Eq. (gni) 

and dehning a) = uj/mb we can write Eq. (El as 


= jdu C^''\ij,fi/mb) , 

= {LW)^K. (53) 

Thus our notation is that contains the Wilson coefficient, while Jq^\uj) is purely the 
held operator. With the convolution notation in Eq. ( 1 ^ the Wilson coefficients are just nu¬ 
merical functions which do not transform under RPI. We will often suppress the dependence 
of Wilson coefficients on ^/mb in what follows. 

Next consider currents that are suppressed by a power of A. At this order the only 
additional structure we can use is a ~ A, where the derivative acts to the left or to the 
right. To form the most general collinear gauge invariant we take D/r^W), which we 
then insert between the {inW) and hy to satisfy the usoft gauge invariance. Since the two 
collinear factors are invariant by themselves they can have arbitrary labels Ui^ 2 - Thus we 
have operators with the structure 

(Lwur T h,. (64) 

The factor of 1/P^ is included to make the Wilson coefficients dimensionless. To make a 
scalar current the fi superscript in Eq. (El can be dotted into a 7 ^ or in P. In either case 
the most general remaining Dirac structure involves either 1 or ^ as follows from section lllBl 
Thus, combining the constraints from gauge invariance, spin reduction, and type-III RPI 
leaves eight 0{X) suppressed currents 

= Jduidu2 b‘f\u;i,u;2) ICj'\uJi,uj2), (55) 
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where j = 1 ,..., 8 , and the bi coefficients are dimensionless functions of 0 ) 1,2 = ^ 1 , 2 /mb and 
^/rrib. The eight operators are 

/cg„Oi,a,2) = 

/egg,Oi, 0 , 2 ) = (inWU{WUv^^WU~[^,^]h,. (66) 

Note that the dependence of the Wilson coefficients on the labels oji account for insertions 
of in all possible locations. Just as for the leading currents we cannot use n-v to form a 
type-III invariant in Eq. dsni) as it leads to currents which can not be made invariant under 
type-II transformations (the transformed currents would depend on e±-v/{n-v)‘^ in a way 
that could not be canceled). 

Next consider the type-I transformations for the currents in Eq. dSSl)- For these subleading 
currents only the A*’ transformations are necessary since we are only working to order A. 
Under type-I only transforms and we have 

= -^iWUn-D,WU = +^VU{u2-n-vV^), 

(WUn -D.WU = -^P 6{uj2-n-vV^). (57) 

Since Vhy = 0 it is easy to see that each of 7 g are type-I invariant all by themselves. 
The other currents do transform, and using Eq. m gives 

jf)/CfHa;i,a;2) = S{u;2) 

6f^ici^\u;uu;2) = bM (enlU)^i(^^)h., 

6f^)C[^\u;,,u;2) = SM > 

*^1 7 , 8 (^ 1 /^ 2 ) = 0. (58) 

The delta functions J(a; 2 ) cause only the coefficients &i *2 3 4 ( 1 ^ 1 ) 0) to appear in the transfor¬ 
mation of A^|i ^2 3 4 }- need the order A variation of the LO current in Eq. In 

this computation we must be careful to note that the 6{uj-n-vP^) in depends on 

n-v, and therefore also transforms 

^ S{uj—n-vP^) = —n-A_i_'P^ S'{LJ — n-vP^) = -- — — 00 6{uj—n-vP^). (59) 

n-v doj 
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Using Eq. (0^ we find a term from transforming the delta function and a term from trans¬ 
forming the collinear quark held 

J‘*V) = A ^K . ( 60 ) 

V 4 / n-v dio 

Demanding invariance under the transformations in Eqs. (EHD and (EOl) gives non-trivial 
constraints on the Wilson coefficients in and From Eq. (EHD the currents 2 3 4 } 

are invariant by themselves provided that 3 4 }('^i) 0) = 0- However, in Eq. (IHUll the hrst 
term can only be canceled by a with 0) = Cq^\u}). To cancel the second term we 

integrate by parts to give a. ujd/duj = ujd/duj acting on Cq \ This term can then be canceled 
by 64 *^ (chi, 0) = —2Cjd/dCj Cq\ijj). Thus, the summary of type-I invariants is Jo-t-ZCi"^^ , 

^{ 5^6 7 8} ’ ^{ 2 ^ 3 } ’ any 65 ®^ 7 g(cDi, 0 ) 2 ) and coefficients 

6i'^^(a),0) = Cq\uj) ^ 0) = 0 ) ^^'’^(ci), 0) = —2(h d/dch (^(^^^(ci)). (61) 


Now consider the type-II transformations. The analog of Eq. (isnD is 

5{uj—n-vV^) = —n-n£_L-P| 5'{uj-n-vV'^) = —n-ve±-V]_-^ 6{u — n-vV^). (62) 

duj 

For the order A variation of the LO current jlf\u!) we have terms from the transformation 
of the delta function, the collinear quark held, and the Wilson line W 


6 


(A) 

II 


Jo{uj) 


duj V in-Dc 2 

hy. 

\ in-Dc 2 u} 


(63) 


In the subleading currents both n^ and have A° transformations.® A ilj)^ trans¬ 
forms to give a so since i/i^n = 0 it is easy to see that ZC| 2 ^g| are invariant under type-II 
transformations at this order. The transformations for the remaining currents are more 
involved 

■SifAg,,(0,1,0,.) = 0, 


4"">k:<”'(o,i.o,J = (LWU(W'i9iW)J-^^K , (64) 

>c<”>(o,i,o,i,) = (ar)„.{A(n/ti„.Spio„ - 
4,U4">(o,i,o,2) = 

4f>4*>(o,i.o,2) = (f»»')o.(irti9p4')„.h^ft„, 

4“>4*>(o,,,o,i,) = [mu{i^{wUv^iWU - ’!:::!^{wHEy:D^WU]~K, 
4">4”>(o,i.o,.) = -\(LWU{WheE^^WU^^K, 


® Since SnV = D(A) the transformation of the delta functions in 



only appear at one higher order. 
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It is straightforward to see that it is not possible to form a type-II invariant using only the 

(s') 

currents 34578 }- However, it is possible to form an invariant taking a combination of 
/C|*i ^4 8 } with Jq\ To facilitate this we rewrite Eq. as 




d 

n-v — 
duo 


1 


i^nWUiW^ (enhE).i 


1 


a;i—a; 


-{WHes_-DjW)^.^,K . ( 65 ) 


To derive Eq. (jnsi) we used e±-V± = ^ni£±- D^Wh^ + ^ni£±- DjoWh^, and the fact 

that the a;i integration can be carried out with the delta function in (^nhE)tJi to get back 
a product of operators with momentum label uo where the intermediate WW^ cancel out. 
Now for the d/duo terms in Eq. (IF)H|l we can integrate by parts in [Co(c5) ^ Jg*^(a;)] so that 
the derivative acts on the Wilson coefficient Cq^\co). It is then evident that the third term 
in Eq. (jnsi) can be canceled by with a) —cDi) = Co*^(a)), the second term is 

canceled by with 64 '^^(a’i, O’—cDi) = —2ujd/dCjCQ\uj), and the hrst and fourth terms 

are canceled by 5^ with 6 g'^^(cDi,ch — cDi) = —2ujd/dcdCjf\Lo) — 2 uo/{loi—uj)Cq\uj). 

(s) 

Therefore, type-II RPI rules out the operators /C 3 5 7 and leaves only the invariants 


+jcr+v:> + jcr , 




{«) 


-(s) 


/c. 


{^) 

2 ) 


/c, 


(^) 


( 66 ) 


For 


type-II invariance their coefficients can have any ^y{idi,uj 2 ), but require 


b[^\ui,u-ui) = C^^\tb ), b[^\u,u-uji) = -2uo -^Cq\u) , (67) 

CLUJ 

b['\u;i,u;-uji) = -2u-^C^'\u;) - . 

CLUJ —UJ ) 

The restrictions on the Wilson coefficients are summarized in Table d Comparing the 
invariants in Eqs. (EQ) and dnzi), we see that the combinations in Eq. are the most 
general combinations invariant under type-I and type-II with the restrictions in Eq. dnzi) 
plus 62 ^^(ci), 0 ) = 0. However, with this constraint on b^^ the operator is actually 
identical to /Cg*^ with an unconstrained coefficient 6 g*\ To see this note that within square 
brackets [WUl^iWU = [WUl^^WU = -[WUg^WU/uo 2 , so the difference comes from 
iljj acting also on {CnW) in . However, since the factor (^nhh) on the left is a collinear 
color singlet we can write 

(w'uftwu = . ( 68 ) 

and the last term vanishes since 62 ( 1 ^ 1 , 0 ) = 0. Given this result and the constraints in 
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RPI-I 

6 f^(+i, 0 ) = 

RPI-II 

6 f^(+l,+ 2 ) = 

Combined Constraints 



Oq ^(+i++ 2 ) 

5 ^^^(+ 1 ,+ 2 ) = C'|^'*^(+l++ 2 ) 

"2 

0 

6 ^®^(+ 1 ,+ 2 ) 

b 2 \uji, 0 ) = 0 

Lb) 

"3 

0 

0 

b'f\ibi,uj 2 ) = 0 

ftS*' 

- 26 )lCi*^'( 6 )l) 

— 2(+1++2) Og ^ (+I++2) 

6^^^(+1,+2) = -2(+i++2) Cjf^ (+1++2) 

i-r 

6|*^(())i,0) 

0 

4''^(+i,+2) = 0 

Lb) 

"6 

(+1,0) 

6[i*^(+l,+2) 

b^Q ^ (+1, +2) unconstrained 

b\^^ 

6^*^ (+1,0) 

0 

b'f\uii,ui 2 ) = 0 

l(s) 

"8 

6^*^ (+1,0) 

—2(+i++2) Cq ^ (+1++2) 

bg ^(+i,+2) = —2 (+i++2) Cq ^ (+1++2) 



T'o*^(+1++2) 

+ S^) C'o^^(+l++2) 


TABLE I: Summary of RPI constraints on the coefficients of the scalar currents in Eq. dSni)- The 
first column shows the constraints from type-I RPI on 0), the second column shows the 

constraint on ( 0 ) 1 , 61 ) 2 ) from type-II RPI and the third column gives the combined constraint. A 

generic entry, like {(1>i,lv2) in the second row of the RPI-II column, indicates no constraint. The 
final currents are displayed in Eq. dzu, and are defined so that they automatically satisfy these 
constraints. 


Eq. (IHTfl it is convenient to define 


= Iduji }Ci\uJi,uj—uJi), 


(s)/ 


K2 {uj) = Cj Jdui {/C4 (coijio;—6i;i) +/Cg (cnijco—coi)} , 

K^^\caJ) = U! JdiOi , 

K^\uji,uj2) = ( 601 + 61 ) 2 )/Ce ^( 1 ^ 1 , a; 2 )- 


(69) 


From RPI it is only these operators that can ever appear. RPI has ruled out some currents 
and restricted R") 2 3 to only depend on one parameter. Once we know this, we can simply 
forget about the and work directly with the Using capital R’s for their Wilson 

coefficients, our hnal basis of subleading scalar operators is 



K'd = 

Jdu K\%(u,), 

\ 


= 

/ <hjiduj2 (+1 , +2 ) {uji,uj2), 





J 


(70) 
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where 



and the RPI type-I and type-II constraints on the Wilson coefficients become 




(s). 


= -2C<‘> V) . = -2C'‘'(ii) 


)(s)/ 


'<(■5)/ 


( 71 ) 


(72) 


The prime here denotes a derivative with respect to a). Thus, we conclude that there are 
4 subleading 0{X) scalar heavy-to-light currents. The coefficient b[^\u}i,u} 2 ) is completely 
unconstrained, while the other coefficients are hxed from RPI invariance. 

In Ref. Q it was noted that Jq\ k[^\ and are connected by RPI, and for these op¬ 
erators our results agree with taking their matching calculation in Eq. (120) and multiplying 
by a common Wilson coefficient. The operator does not appear in Ref. , because 
the derivative on its coefficient causes it to vanish at tree level. Our three-body operator 
is also new. In the limit that depends only on the sum cD+ = a>i -t- 0)2 we 

can switch variables to cD+,ci)i and reduce to a two-body operator. At tree level this 
is always possible since the Wilson coefficient is independent of the ooi. To see how the 
reduction works we write 

[du,Ki^\u;i,u;+-uJi) = —{LW^W^ig^w) K 

J rUb '• V /uj+ 

^ —{L[i$iW]) K, (73) 

rUb \ / 07 + 


where in the last line we used Eq. flldp. T he derivative structure of this two-body operator 
is similar to that of operators in Ref. [l^ ll^ , however the specihc spin structure appearing 
in Eq. dZSl) does not appear from matching the QCD scalar current at tree level. Beyond 
tree level R 4 can depend separately on uji and a ;2 and the reduction in Eq. dZSD is not valid. 


B. Vector Currents 

The steps for deriving the general set of vector currents are very similar to the steps for 
the scalar currents in the previous section, so our presentation will be more concise. 

At LO gauge invariance plus power counting allows currents of the form h^. Im- 
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posing type-III RPI invariance allows the Dirac structures 


Z 




rinU-v, —^ 

n-v 


n-v 


, n^^ , 


n-v 


n^^{n-vY 


n^yfin-v 

n-v 


\ 

{n-vY / ’ 



n-v ’ 


( 74 ) 


However it is easy to show that all structures involving n are ruled out by the % trans¬ 
formations. Thus, at leading order there are only three allowed vector currents 


T-h) 

"'l-S 


chi c!”_> 3 (i) 


(v) 


(75) 


where the coefficients are functions of cD = a;/mb and /i/ruft, and 

( 76 ) 

The choices {1,2,3} correspond to the three different Dirac structures, and our basis in 
Eq. (|7F)|l agrees with Ref. 

At NLO the power counting only allows a single D± to appear. For the possible spin 
structures it is easy to see that type-II RPI invariance does not allow the vector index to 
be in an n^ or any factors of n-v to appear just as for the leading currents. Imposing the 
constraints from gauge invariance, spin reduction, and type-III RPI then leaves 28 0{X) 
suppressed currents [j = 1 ,..., 28] 


ICf'^ = ldujidLU2 h'f’{id;i,uj2) ^^ 2 ), 


' h) 


h), 


( 77 ) 


with operators 


/Ci%(o.i,a;2) = . 


/c!l 2 o(^i,^ 2 ) = {Lwu{w^iij)^wu 


pt I 2 ’n-n 
1 

pt 

1 


hv 7 




~ — \K, 




pt 12 n-v 




Pt 


where the six P(^ matrices are 


■p/7 f ^7^ 7^ 

^ 1-6 “ 1 O ’ o ’ ’ 


n^ 


I 2 ’ 2 ’ 2n-v ’ n-v ’ n-v ’ {n-vY 


( 78 ) 


( 79 ) 


Working out the transformations of the leading and subleading currents in a similar way 
as was done for the scalar currents we hnd that the type-I invariants are -|- + k!Y 2 1 
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RPI-I 

= 

RPI-II 

6f4‘ai,ch2) = 

Combined Constraints 

'^1-3 


C[\{cbl+UJ 2 ) 

6i%(o;i,(h2) = C'{%((hi-hch2) 

ffi-6 

0 

4-6(‘^1i‘^2) 

626(^^11 0 ) = 0 

Oj 

0 

0 

b 4 iibl,ib 2 ) = 0 

b^) 

- 2 c 4 {C,r) 


6f((hi,0) = -2Cf4(hi) 

4-11 

0 

0 

69!^ll((hi,(h2) = 0 

'^12,13 

—2(hi i'^i) 

— 2((hi-1-0)2) (o)i-|-(h2) 

61273 (“^li *^2) = —2(10)1-1-0)2) (161+^2) 

"14 

- 2 ( 01 ^ 4 ' 

—2((hi-|-ch2) C*3 ^ (chi-|-(h2) 

644"^1i"1’2) = — 2 {ll)i+ill 2 ) C3 ^ (chi-1-0)2) 



-6^'’4wi,h)2) 

-6^'’4wi,(i)2) 

"15-17 

^15-17(^1’ b) 

0 

64-17(^1 i "I’s) ~ b 

"18-20 

^18-20(^^11 b) 

^18-20 (^li ^2) 

^18-20(‘I"!’ ‘^2) unconstrained 

"21 

4?(‘^iib) 

0 

6^3^ ((hi, (02) = 0 

AA 

"22 

^22 (‘^ij b) 

4?(‘^1>‘^2) 

b^22 unconstrained 

"23-25 

43-25 (‘^li b) 

0 

64-25(^^11^2) = 0 

A^) 

"26,27 

4677(^^11 b) 

—2((hi-|-ch2) C*jJ (o)i-|-(h2) 

6477(^11^1^2) = —2((hi-1-0)2) (71^2 (161+^1)2) 



—2 (7^2 (‘-"1 +^2) 

—2 C^2 (^1 +0)2) 

00 

-0 

4? (‘^lib) 

—2((hi-|-ch2) C*3 ^ (chi-|-(h2) 

644"^1i" 1’2) = —2((hi-|-ch2) Cg ^ (chi-1-0)2) 



+ Sj) C'3^4‘^1+‘^2) 

+2(f + Sj) C'3'’4o)i-h(h2) 



— 622(0)1,(1)2) 

—622 (0)1,162) 


TABLE II: Summary of RPI constraints on the coefficients of the vector currents in Eq. dZHl). 
The first column shows the constraints from type-I RPI on b^^\LO, 0 ), the second column shows the 

constraint on b\^\L0i,uj2) from type-II RPI and the third colnmn gives the combined constraint. 
The final cnrrents are displayed in Eq. m, and are defined so that they antomatically satisfy 
these constraints. 


-^2 '' + ^2 '* + ^13^) ^ + ^3 '* + ^8 + ^14 ) ^'j ’ ^ ^ where j = {4, 5, 6, 7, 9,10,11} 

and k = {15, ...,28}. Type-I invariance allows any coefficients bnt restricts 

^1-14 ('^1;'^2) as shown in the second colnmn of Table ITTl 

Looking at the transformations nnder type-II we hnd that the invariants are independent 
of the {7^, Vfj_, n^} choice. Onr resnlts for the type-II invariants are , 


-(v) 


-h) 


-h) r-h) 


(v) 


4 '"^ +164 + ^ 

-h) 


h) 

13 


+ /C 


h) rh) 

27 ! "'3 


+ /C<”> + /Cf' + /C” + (K™ + K'S ), (K™ - /C£’), while 


-h) 


-h) 




(v)\ 




-h)^ 


KJf for £ = (4, 5, 6,18,19, 20} are invariant by themselves. For these combinations type-II 
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invariance allows any 0 ) 2 ), but restricts 9 - 17 , 21-28 shown in the third column 

of Table |TT1 Furthermore, currents Km with m = {7, 9,10,11,15,16,17, 21, 23, 24, 25} are 
ruled out (ie. bm'^ = 0). 

It is easy to see that the type-I and type-II conditions in Table HH are compatible. The 
combined set of constraints are given by those in the fourth column. Using Eq. dHED we 
can show that the constrained K^^^^ g} are redundant with /Cj^g 29 } respectively, just as 
was done for the scalar current with K^'^ and Kq\ We can also use Eq. (IF)H|) to convert 
K^\( 2 !i,u} 2 ) — 7 F| 4 ^(a)i, 0 ) 2 ) into a term proportional to 6 ( 002 ) and a term that is redundant 
with K^I\o 6 i,l 62 ) - K^s ( 661 , 002 )■ 

From the combined constraints we can then dehne a new complete set of allowed vector 
operators, Therefore after imposing type-I and type-II RPI plus all other constraints 

we are left with our hnal set of allowed vector current operators 



_ 

-f^i-io ~ 

jdu, 




-'^11-14 ~ 

Jdooidoo2 RnLi4(a)i,a;2) iFn)_i4(a;i, 012 ), 


(80) 

\ 



J 



where 



(81) 


The coefficients 42 , 13 , 14 } ^9- dSOl) depend on two parameters 011^2 and are uncon¬ 

strained. The remaining coefficients depend on only one parameter and are fixed by repa¬ 
rameterization invariance 


= cr4(£.), Br-7(^) = -2cr4'(£.), = - 20 ^ 0 ,^ , 




){U 




)h) 






h)/ 


(82) 


The form of the currents agree with Ref. 3|, and if we take a frame where n-v = 1 
and = 0 then also agree. Ref. US looked at type-I RPI of the vector currents 
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and our constraints on agree with the ones found there. (We note that the authors 
of Ref. l3| also checked these results with explicit one-loop computations.) At tree-level 
one matches onto the currents K^l and the two-body limit of (using the analog of 
Eq. ((ZSI)), and we agree with Ref. on the form of these currents and the RPI constraint 


between and The structures in Eq. ea which are new and which only appear 


.{v) 


beyond tree level are and the three-body form of . 


(v) 


C. Pseudoscalar and Axial-Vector Currents 

The results for the pseudoscalar and axial-vector heavy-to-light currents can be directly 
obtained from the analysis for the scalar and vector currents respectively. The analysis is 
identical except for the extra 75 in the Dirac structure. For the pseudo-scalar currents we 
have the basis {^/ 2 , l/n-n} 75 , while for the axial-vector currents we have r|) g| 75 where 
Tj is defined in Eq. dzni). At LO the most general allowed pseudoscalar current are thus 

= jdu Co^\u,fi/mb) Jo^\u), 

(83) 

while the axial-vector currents are 

J<7 = jdu , 

jfcUo;) = (84) 

At NLO we again have eight possible pseudo-scalar currents K.^ and 28 possible axial- 
vector currents K.^ before imposing all type-I and type-II constraints. After imposing the 
RPI constraints we find results very similar to those in Eq. (ED) and m- Thus for the hnal 
NLO pseudoscalar currents we have 



Aty = jduj 

CO 




= [duiduj2 

(<hi,ch2) A"4^^(a;i,a;2) , 


( 85 ) 

V 



) 



where 






and the RPI type-I and type-II constraints on the Wilson coefficients are 

= -2C‘^\u). (87) 

For the hnal axial-vector NLO cnrrents we hnd 








-'^11-14 ~ 

Jdujiduj2 Sn^_i4(a)i,ch2) iFnl_i4(a;i, a;2), 


(88) 




J 



where 



The coefficients 13 14} Eq. (jHn|l depend on two parameters a;i,2 and are nncon- 

strained. The remaining coefficients depend on only one parameter and are hxed by repa¬ 
rameterization invariance 

= -2C^^\Cj). (90) 
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The form of the pseudo-scalar and axial-vector currents are very analogous to the scalar 
and vector currents, so rather than comparing with the literature we simply refer to the 
comparisons in the proceeding sections for which part of our results were previously known. 


D. Tensor Currents 


At leading order in A, there are four tensor currents, defined as 

= {(nWl , 

where the most general allowed Dirac structures are 






(91) 

( 92 ) 

n-v n-v 

where = 7 ^nj, — 'juVfj, etc. As before, no can appear at leading order from type-II 

RPI. 

At 0(A), 44 currents can be written down before imposing the RPI constraints. They 
can be chosen as 




V 


(93) 


4-14 = (LwUiwHDl^^wUTt,^ K , 


and 


/C 

K, 

/C 

/C 


V 

S -22 = (UWU (ir* iv^^wu \ 4^ a. , 

^ iL’U j 'jj? 


(b 

23-30 


= (f„tr)„,(W't i$iwu I hti . 4-rr41A K , 


n-v 


V' 


(94) 




31-36 

(b 

37-44 


V 


= {^nWUiw^ tv-Djwu I|rr 4 , — rr 4 l i 


2 'n-v 

The Dirac matrix with one index T'^_q is defined as in dZHD 

1 / n,, 


V' 


^1-6 = { I — 

' 2 \ v-n/ n-v 


7/4) 


v-n 


(96) 


The constraints from type-I and type-II RPI are derived as before. The final constraints on 
the Wilson coefficients 61-44 are shown in Table uni 

After imposing the constraints from the table one finds the final minimal set of tensor 
heavy-light currents in the effective theory at 0(A) 


( 96 ) 




N 

^(b _ 

-'^1-14 ~ 

1 da,B!7,(^)/f!‘l„(^), 


^(b _ 

-'^15-21 ~ 

da;idu;2-S|5L2i('^i) '^2)A"f5A2i(i^i) 1 ^ 2 ). 


V 
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RPI-I 

6 f^(chi, 0 ) = 

RPI-II 

bf\cbi,cd2) = 

Combined Constraints 

^(i) 


C[}_^{cbl+ib2) 

45.4 (± 1 , ± 2 ) = C'f 24 ((hi+(h 2 ) 

^(i) 

°5-8 

0 

bfl^{Lbi,Lb2) 

4 l 8 ((hi, 0 ) = 0 

bit) 

"9-11 

0 

0 

45.11 (± 1 , ± 2 ) = 0 

bit) 

"12 

2 Cf((hi) 

65 §((hi,(h 2 ) 

45 (± 1 , 0 ) = 2 61^ ((hi) 

bit) 

"13 

- 2 ^^ ((hi) 

^13 (‘^ 1 ;‘^ 2 ) 

6 g((hi, 0 ) = - 2 Cf(±i) 

bit) 

"14 

0 

6 f]((hi,(h 2 ) 

45 (± 1 , 0 ) = 0 

1 

00 

0 

0 

45-18(^1 >± 2 ) = 0 

^19,20 

—2 (hi cj 2 {‘^ 1 ) 

— 2 ((hi+(h 2 ) C *^2 (^lT^ 2 ) 

^194o(±1’±2) = — 2 ((hi+(h 2 ) (±l+± 2 ) 

bit) 

" 21,22 

—2 (hi C 3 4 ((hi) 

— 2 ((hl+(h 2 ) ((hl+(h 2 ) 

4 ^ 22 (± 1 ’± 2 ) = — 2 ((hi+(h 2 ) ((hi+(h 2 ) 


+2 C'3!4(‘2’i) 

±45,13 

±^ 1243 (^ 1 ’ ± 2 ) 

bit) 

"23-26 

^23-26 (^1’ 

0 

^4426 (± 1 > ± 2 ) = 0 

^27-30 

^^-30(^1’^) 

47-30 (± 1 ’ ± 2 ) 

627 _ 3 o (± 1 , ± 2 ) unconstrained 

4i-33 

4 i- 33 (‘ 2 '!, 0 ) 

0 

4i433(±1,±2) = 0 

^34-36 

^4-36(‘2'!, 0 ) 

44-36 (± 1 )± 2 ) 

45-36 (±1, ± 2 ) unconstrained 

^37-40 

^^-40(^1’^) 

0 

4744o(±1,±2) = 0 

bit) 

"41,42 

^41,42(^1) 0) 

—2((hi+(h2) (±lT±2) 

4p42(±l’‘^2) = —2((hi+(h2) (±l+±2) 



CiS(±i+±2) 

+2 (1 + ^) C'p2(±1+±2) 

bit) 

"43,44 

^43,44 (^1 5 

— 2((hl+(h2) ((hi+(h2) 

43,44(±1’^^ 2 ) = —2((hi+(h2) (±3^4 {UJ 1 +QJ 2 ) 



+^(1 + ^) C'44(±1+±2) 

±'3!4(‘^1+‘^2) 



±^34,35 (h>l,±2) 

±^34,35(±1, ± 2 ) 


TABLE III: Summary of RPI constraints on the coefficients of the tensor currents in Eq. (P|). 
The first column shows the constraints from type-I RPI on bi{to, 0 ), the second column shows the 
constraint on 6j(wi, CU2) from type-II RPI and the third column gives the combined constraint. Each 
± refers to the first and second terms in their row respectively. The final currents are displayed in 
Eqs. (Ezm, and are defined so that they automatically satisfy these constraints. 
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There are 14 independent 0(A) two-body operators given explicitly by 



Their coefficients are hxed by reparameterization invariance in terms of the 0(A°) Wilson 
coefficients as 

£?|‘y(w) = c|'y(w), = 2Cf (OJ) , BfV) = -2Cf(w) . 

= -2C[‘-4(‘^) ■ (98) 

In addition, there are 7 three-body collinear operators given by 



and have to be determined by an explicit matching calculation. 

At tree-level one matches onto the currents j and the two-body limit of (using 
the analog of Eq. CHI), and we agree with Ref. jl^ on the form of these currents and the 
RPI constraint between and The remaining operators in Eqs. (IHZI) and (jnni) are 

new and only appear beyond tree level (including the three body structure of ). 


V. SUMMARY FOR COEFFICIENTS, OPERATORS, AND FEYNMAN RULES 

In this section we summarize results that should be useful for future phenomelogical 
applications. In Section IV A I we summarize the full set of known matching results and 
compare with the literature, in section IV Bl we give simplified expressions for our basis of 
currents in the frame = 0, n-n = 1, and in section ED we give Feynman rules for the 
subleading currents and ■ 
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A. Matching results for the currents 


In this subsection we summarize the one-loop matching results for the LO and NLO Jhi 
Wilson coefficients. With = uj/rrii, these coefficients are defined in previous sections as 
and ^/rrih) respectively, where (d) denotes whether the current is a 

scalar, pseudoscalar, vector, axial-vector, or a tensor. 

For the LO currents the basis we use is different (though equivalent) to the basis used in 
Ref. . Since the one-loop matching for the LO coefficients can be found in Ref. [3 it is 
useful to have the explicit relation between our basis of coefficients and the coefficients 
Cj (j = 1-12) that can be found there. We find 


'^0 

= c,, 

C^v) 

= 03 , 

OjO 

= 06, 

of) 

Mp) 

'^0 

= 02, 


= 05 , 

of) 

= 08, 

of) 




= 04 - 03 , 

of) 

= 07 - 06, 

of) 


of 


Oio, (100) 

—012 , 

Olo-Og, 

O12+O10 —Oil • 


At tree level the matching between QCD and SOFT is scheme independent. Matching with 
the full QCD currents -u{l, 75 , 7 ^, 7 ^ 75 , we find 


_ r''-p> — n''’") — — 1 

'-'0 — '-'0 — '-'1 — '-'1 — '-'1 — ) 


dp) 




(a) 




_ ^ya) _ ^yu) _ ^yu) _ ^yLj _ ^yi) _ „ 

0^2 — *-^3 — '-'2 — *-^3 — '-'2 — '-'3 — '-'4 “ • 


'<0) 










di) 


( 101 ) 


At one-loop we use the MS scheme with naive dimensional regularization (NDR) and match 
at fi = rrib to determine the fi/mb). Using Eq. (Il()()|l and results in Ref. [3 one hnds 


O, 


is,p)i 


0 ( 0 , 1 ) = 1 


ofO(^^l) = 1_ 

of (0,1) = 1- 

Of0(^,l) = 


as{mb)CF 
dvr 

as{mb)CF 
Air 

as{mb)CF 
dvr 

as{mb)CF 


21 n^(a>) + 2 Li 2 (l— 0 ) 


2 ln(a)) TT^ 
1-0 12 


21 n^( 0 ) + 2 Li 2 (l— 0 ) - 1 - ln( 0 ) 
21 n^( 0 ) + 2 Li 2 (l— 0 ) -f ln( 0 ) 


30-2 

1-0 

dO-2 

1-0 


71 




TT 

H-h 6 

12 


dvr 


+ 


20 ln( 0 ) 


(i-i,) (i-i.)q• 


7 ‘>(ii,i) = 0, 


c7“>(ii,i) = 
Cf(i,l) = 


Q;s(m;,)Oir f (1 — 20)0 ln( 0 ) 0 


dvr ( (1 — 0)^ 

as{mb)CF f-201n(0) 


u 


dvr 


u 


Cf(ii,l) = 0, 


( 102 ) 


where Cf = d/3 for color SU{3). To determine the coefficients for scales mfoAgco < < ml 

we require their anomalous dimensions.® The LO and NLO anomalous dimensions are 


The full NLO result requires a two-loop anomalous dimension which uses information from Ref. |3 !t|1 
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universal and the running of these coefficients is given in Ref. (or for the case cD = 1 in 
Ref. Q). 

At NLO in A tree level matching of the QCD current uTh onto SCET gives TO, H, 3 
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7(16) 

^tree 


7 ( 1 ^) 

^tree 


= -Wr| 

0 _ 1 


ho) 


n-v 


in- D, 


n-v mb 


iv-D^WK 


(103) 


The normalization in Eq. (unsi) was hrst derived for in Ref. H. and for in 

Ref. jl^. Comparing Eq. (Unsi) with our basis of currents we see that for any choice of E 
the Jtree^^'^ match onto a subset of our two-body currents. On the other hand J^^el does not 
appear in the basis of two-body currents. Instead it is obtained from the projection of a 
subset of the three-body currents for cases where the corresponding coefficients Rj( 0 ) 2 , 1 ^ 2 ) 
depend only on the sum a)i -|-ci> 2 . This is certainly the case at tree level since the coefficients 
are Ui independent. The three-body structure of the currents can only show up at the level 
of one-loop matching. 

Using Eq. (uni to determine the tree-level value of the NLO Wilson coefficients of the 
operators in Eqs. (ITWn El EHl EH EHI) we End 
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-°14 
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-°ii 
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oh) _ 
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(104) 

These results are in agreement with the RPI constraints in Eqs. dZl El Ell El El)- Coeffi¬ 
cients in Eq. dnn that are zero indicate that the corresponding currents can not be inferred 
at tree level since they are hrst matched onto at one-loop (or beyond). The full one-loop 
matching for all the 0{\) currents is not currently known from direct computations. How¬ 
ever, many of the NLO coefficients are hxed in terms of the LO coefficients by RPI, namely 
Ro- 3 ^(<^)) and Rf2]^4(a>). Summarizing Eqs. ()72l82l87l9()l98jl we have 
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11-14 
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(106) 


where to save space the cD and /r dependence of the expressions on both sides of these equali¬ 
ties is suppressed. These results can be used to determine the matching for these coefficients 
at fi = mb using Eq. (mn). They also imply that the anomalous dimensions of these coef- 
hcients are determined by the anomalous dimension of the leading order coefficients j^, so 
their values for scales mfcAqcD < < mf is known. 
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For the coefficients of the 3-body operators, ^ 2 ), and 

i?|; 5 _ 2 i(cDi, 0 ) 2 ), neither the one-loop matching results, nor even the LO anomalous dimensions 
are currently known. 

Finally, we note that it is possible to relate the pseudoscalar and axial-vector coefficients 
from the scalar and vector coefficients. For massless quarks the QCD diagrams and SCET 
diagrams change in a trivial way under the chiral transformation, q —>■ 75 ^, and 
provided we work in a scheme such as NDR. Therefore in this scheme the Wilson coefficients 
of operators with and without 75 are related (see for example Ref. for the relations between 
LO coefficients). In other renormalization schemes these coefficients may differ. 


B. Summary of 0{X) currents in the frame v± = 0, n-v = 1 


In sections IIV Al through IIVDI we have derived the most general basis of heavy-to-Iight 
current to 0{X) in an arbitrary frame. However, for applications it is often most convenient 
to pick a frame where n_L = 0 and v-n = 1. In this frame the currents K 2 f \ and 

drop out. Thus there are only (2,2,8,8,13) order 0{X) heavy-to-light currents which 
are (scalar, pseudo-scalar, vector, axial-vector, tensor). In this section we summarize our 
results with this choice of basis vectors. 

In this frame our leading order results for the Jhi currents with a complete set of Dirac 
structures can be summarized as 


j 7 ( w ) = i.UW)^TfK, ( 106 ) 

where in F-^^ the (d) specihes the type of current (scalar,vector,...), the i specihes the 
member of the complete set of possible structures of that type, and the Wilson coefficients 
are For the minimum basis of Dirac structures we found 

p(®) 1 p(p) 

-*-0 to ~ 75 5 

r{i!2,3} = {7 m 
^{ 1 , 2 , 3 } {TaiTs ) 

^{ 17 , 3 , 4 } = 1 d[M,Wl 1 d[M,W] , , (107) 

which is simply a linear combination of the basis in Ref. [^. At order A our corresponding 
results for Jhi in this frame can be summarized as 


j(la) _ 


(108) 
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(f„ iSy w)^ d t 7» /i„ , 
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* n-q 


{d)a - 


FIG. 1: Feynman rules for the 0(A) currents in Eq. (nnni) with zero and one gluon (the fermion 
spinors are suppressed). For the collinear particles we show their (label,residual) momenta, where 
label momenta are p, ~ and residual momenta are k,t ~ A^. Momenta with a hat are 
normalized to mf,, p = p/rrib etc. 
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■2 
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.275, 

0 (P)« , 

'-^11-14 ~ 

{7"72, 
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-' 11-14 ~ 1 

7''757^ 


(109) 




qW« _ 


0 


[fi u] a [fi iy] a a[fi v] a[^i v] a[^i u] 


15-21 = i 7 ± , T’'v '1± , r^'n ^7^ , gp'-'q ' , gp'-’v ' , gp'-'n 


Note that due to Eq. (HI the form of in Eq. dniHi) is identical to the form of the 
currents that was used in Ref. 0 since [hF'^i D^^W\ = [l/P WHgB^^W]. 


C. Feynman rules for J^i and C^q 

In this subsection Feynman rules are given for the 0{X) heavy-to-light currents Jd“) and 

J( 16 ) 

in Eq. (HUHi) which are valid in a frame where u_l = 0 and v-n = 1. We also give the 
Feynman rules that follow from the hnal form of the £d> 2 a, 26 ) L^grangians in Eq. dH. 

For the subleading currents the zero and one gluon Feynman rules for Jd“) and are 

shown in Figs. ^ and |21 respectively. (From the results in the previous sections the Feynman 
rules for the currents with u_l 7 ^ 0 and v-n ^ 1 can also be easily derived.) For the 

Wilson coefficients depend only on the total A° collinear momentum, while for the 

coefficients depend on how the momentum is divided between the quark and gluons. The 
j(ia) cai-rent has non-vanishing Feynman rules with zero or one gluon and any number 
of n-An gluons. The possible gluons that appear in the currents are similar, but the 
current vanishes unless it has one or more collinear gluons present. 
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p, a 

FIG. 2 : Feynman rules for the 0 (A) currents in Eq. with zero and one gluon. For the 

collinear particles we show their (label,residual) momenta, where label momenta are p, q, Qi r\j AOd 
and residual momenta are k,t ~ A^. Momenta with a hat are normalized to m^, p = p/mi, etc. 
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FIG. 3 : Feynman rules for the subleading usoft-collinear Lagrangian with one and two collinear 
gluons (springs with lines through them). The solid lines are usoft quarks while dashed lines are 
collinear quarks. For the collinear particles we show their (label,residual) momenta. (The fermion 
spinors are suppressed.) 


For the mixed usoft-collinear Lagrangians from Eq. dH, 

4? = + h.c. , 

in-Dc 

= In ig]^ W qus + h.c., 

in-Dc 

^ 

all Feynman rules involve at least one collinear gluon. From we obtain Feynman rules 
with zero or one gluons and any number olfi-An gluons. The one and two-gluon results 
are shown in Fig. IHl 

For we have Feynman rules with zero or one {n-An, ^us} gluon and any number of 
n-An gluons. The one and two-gluon results are shown in Fig. Finally, for one hnds 
Feynman rules with zero, one, or two A:j^ gluons and any number of n-An gluons. In this 
case the one and two gluon Feynman rules are shown in Fig. El Note that it is important to 
treat the contributions from and separately since they show up in different parts 
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FIG. 4: Feynman rules for the O(A^) usoft-collinear Lagrangian with one and two glu¬ 

ons. The spring without a line through it is an usoft gluon. For the collinear particles we show 
their (label,residual) momenta, where label momenta are p, q, qi ~ and residual momenta are 
A:, t, ti ~ A^. 
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+n^,n^ ( 

' [n-pY n-pn-q2 


n-q 2 
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FIG. 5: Feynman rules for the O(A^) usoft-collinear Lagrangian with one and two gluons. For 
the collinear particles we show their (label,residual) momenta, where label momenta are p, g, qi ~ 
A'^A and residual momenta are k,tYi ~ A^. 


of the heavy-to-light factorization formulae derived in Ref. m and shown in Eq. ([Q). 

For -|- the Feynman rules are different than one would derive using the in¬ 
termediate form Eq. dSOl), since in transforming to the hnal form the equations of motion 
were applied. However, observable predictions that are consistently made with either set of 
Feynman rules will agree. 
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VI. LEADING POWER PREDICTIONS FOR B TO PSEUDOSCALAR MESONS. 


As a phenomenological example, we consider the form factors for B —> niu, or more 
generally the form factors for B ^ P where P is a pseudoscalar meson (calculations using 
the factorization theorem in Eq. m for the vector meson are just as straightforward). For 
pseudoscalars there are three form factors in QCD, which are conventionally defined by 




Pb+P - 

2 


- mp a 

y 






{P{p)\qia^^''q^h\B{pb)) = - ^ [q'^{Pb+lP) “ ("^1 “ "^p) , 

vriB + mp 


( 111 ) 


A2 

^'■QCD; 

^qcd)- 


where q = pb — p- 

For the region where Q = {E, m^} S> Aqcd (ie. small g^) one can use large energy 
factorization to study the form factors. For pions our expansion parameter Agcn/h-p ~ 
0.5 GeV/(2P) becomes 1/4 for P ~ 1 GeV. This makes the region of where the expansion 
is valid roughly 0 < g^ < 10 GeV^. In SGET the form factors /+, /o, fp split themselves up 
into contributions associated with three momentum regions: Wilson coefficients for 
two jet functions Ja^b for p^ ~ QAqcd, universal light-cone wavefunctions for p^ ^ 
and a single non-factorizable form factor C,p (containing both p^ rsj QAqcd and rsj 
The leading contributions therefore split into factorizable (F) and non-factorizable (NF) 
contributions. This decomposition was dehned by the proof of a factorization formula for 
these form factors in Ref. 

;ecD(^2) _ ^ ^ ^ ^ 

pi pi POO 

f^{Q) = No dz dx dr+T{z,Q,po) 

Jo Jo Jo 

xJ(z,x,r+,Q,iia,li) 4>p{x,ii) 4>B(r+,p , 


( 112 ) 


where P{Q) Q and the ellipses denote terms that are suppressed by more 

powers of l/Q. Here 05 = 0|- 

To separate the scales and QA we match QGD onto an SGETj. The scales QA and A^ 
are then separated by matching SGETi onto an SGETn Q. Operators in SGETi are divided 
into F and NF categories depending on the form of the result of factoring usoft gluons from 
collinear helds. In SGETi the F contributions are from the time-ordered products 


T[ = /‘d^a;T|j('“)(0),2£«(a;)}, 


T4 = (TxT 


1 iPfqix) \ , 


T( = d^xT 


1 iPfa\x) \ , 


(113) 


Here the phrase non-factorizable simply refers to the fact that the matrix elements can not be expressed 
in terms of convolutions with the standard light-cone wavefunctions. 
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where the currents are taken from Eq. (uniD and the usoft-collinear Lagrangians from 
Eq. (j1 lOjl . After factorization of usoft and collinear helds, the Tf^s are matched onto soft- 
collinear SCETn operators. The collinear matrix elements are given in terms of jet functions 
J{z, X, r+, Q) , and the soft operators are given in terms of B light-cone wave functions (f)±{r+) 
dehned as 1110 


{0\q{x )Sn{x ,0)Th^{0)\B{v)) 


(114) 




dr+ 


Tr 


l+j 




75 r 


A few general properties of the factorizable term f^{Q) can be given without an explicit 
computation. First, the matrix elements of T /2 can only depend on 0^(r+). This follows 
from the explicit form of the subleading Lagrangian , where so the usoft held 

Qus appears only in the combination qus^- Using (I114j) this implies that only the 0j(r+) term 
gives a nonvanishing contribution. On the other hand, the factorizable operator depends 
on the combination so its matrix element can only contain 0jj(r+). At tree-level the 
jet function J from the matrix element of vanishes, but a nonzero result could appear 
at one-loop order. However, the matrix element of contains the leading order current 
JW, so it obeys the same symmetry relations as those derived for the nonfactorizable part 
/""«) 0 . Therefore, although this matrix element is factorizable it does not increase the 
number of unknown non-perturbative functions since for phenomenological analyses 0^ can 
be absorbed in . With this choice, all remaining factorizable contributions are expressible 
in terms of just 0^(r+). 

Using the approach explained in 0 we can obtain the results for the form factors. After 
factorization of usoft and collinear helds the T-products of collinear helds coming from Tf 2 
are given by (using Eqs. p()8j) and (Umi)) 
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UJl,UJ2 
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ijB 
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(115) 



iA 

- !jJ2 

( 0 ) 




iJB 

. ( 

J 0 


where i, j are Dirac indices and here A, B are color indices in the fundamental representation. 
The functions J7ia,ife are collinear gauge invariant and satisfy the spin structure constraints 
iJia,ib = Jia^ibi = 0, and ii[Jia,ib] = 0. 

Taking into account constraints from the Dirac structure of the ehective theory helds one 
can easily hnd the most general form of the operators appearing in the matching of 
onto operators in SCETn. The jet functions Ja,b are dehned by the terms which contribute 
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on a pseudoscalar state 




J^^{x) = i5{x+)5 (xx) - Mv . 2 ^ 


e-^r^x 


(116) 
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X Jb(a;,?7,r+) 


(Lw)m-vAy(w'u) 


+ •■■ 


II 


where a; = a;i —a; 2 . We have suppressed the dependence of Jia,ib on the /i’s and on a;, a;i +a ;2 
(the latter combinations would be simply set to n-p in the pseudoscalar matrix element by 
momentum conservation j^). The ellipses in Eq. (HinD denote color octet terms and other 
operators which do not contribute for a pseudoscalar meson P. 

Using Eq. (HinD the operators in {Pn{j))\T[ 2 \Bv) factor into a product of matrix elements 
that can be evaluated with Eq. (imi and Eq.(12,13) of Ref. 0. Switching variables to x, z 
by using uj = { 2 x — l)n-p and p = { 2 z — l)n-p we hnd the following factorization theorems 
which are valid at leading order^^ in 1/Q and all orders in as 
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it). 


where E = E/mb, Q = {E,mb}, and the normalization coefficient is given by No = 
fpfpmB/idEy. The matrix element involving non-factorizable operators gives C^{Q,p) 


11 


We kept a kinematic factor of mp in the prefactor of fp even though it is formally power suppressed. 
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which is the reduced form factor describing decays to a pseudoscalar meson P. The quan¬ 
tities in square brackets and curly brackets are calculable, that is the Ta^s and Ja,b^ have 
expansions in as{Q) and as{\/Qh) respectively. Note that the Ja,b are universal, meaning 
that at any order in as it is these same jet functions which appear for any pseudoscalar 
meson and independent of which form factor /+,o,r we consider. Therefore, the factoriza¬ 
tion theorem still gives information even in the case where we assume that is 

non-perturbative. 

Working at O{as{fio)) (ie. tree level) for the jet functions gives 


Jia{x,r^) 


Jib{z,x,r^) 


nCp (Us (ho) 

Nc a;r+ 

7rCp Oisi^fJ^o) r/ \ 


(119) 

( 120 ) 


At this level the z integrals in Eq. (HIED disappear because the tree level jet gives a 6{z — x), 
and this causes the z variable in the T^’s to be replaced by x. The are combinations 
of Wilson coefficients appearing in the currents given in Eq. (uniD and should be 

evaluated at a scale /ig ~ QA. Expressed in terms of the Wilson coefficients defined in 
Sec. IV, they are given by 


r<+>(J!.i.) 

rf'(B.A) 


B|’>(2£:,ri + 

b|”>(2£:,,^) + 
b!‘>( 2 E,m)- 


(2 B - tub) 

{rriB - 2,E) 

B^•\2E,^L) - 2Bf{2E,p, + B<‘>(2B,;,), 


( 121 ) 


where the dependence on = 2E is shown, and 


tI°'’(e,z,^) 

T^(E,Z,^L) 


B^i>{2E,z,,z) - —Bj”>(2B,i,ri - Bj”>(2£), z,,^), (122) 

tub 

(2E, z, p) - ) (2E, z, /i) - i?;? (2E, z, /i), 

tub 

Bil (2E, z, /i) + (2E, z, p) - sg (2E, z, p), 


where Cji + Cj 2 = 2E and the dependence on z is induced from the cfii — a >2 dependence of 
the coefficients. 

If we work at tree level in Ja^ using Eq. (Illbj) and also in Ta^ then these coefficients are 
scale independent and satisfy = 1 and =0. In this case if we take 

the ratios /o//+ and fp/fp and expand assuming that the terms are smaller than the 
terms then our results agree with Ref. H. We note that using just the information 
in our factorization theorem that it is not clear that one wants to expand in this way since 
the F and NF terms could actually be similar in size as discussed in the introduction. The 
expectation from QCD sum rules is that the “soft” NF part of the form factors is larger 
than the “hard” F part 
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VII. CONCLUSION 


The soft-collinear effective theory (SCET) allows a rich structure of allowed operators at 
higher orders in the expansion parameter A. In contrast to simpler effective theories, the 
presence of fields {n ■ An) and derivatives (h • iDc) scaling like A°, allows a continuum set 
of operators at any given order in A. A similar situation is encountered in deep inelastic 
scattering, where an infinite number of operators of increasing dimension can contribute to 
the same order in 1/Q. In a generic process with energetic hadrons it is therefore important 
to have a well-defined procedure for organizing the structure of the soft-collinear operators 
at a given order in A. This organization is provided by SCET. 

In this paper we formulated a general prescription for constructing the most general 
ultrasoft-collinear operators appearing in the Lagrangian or in the matching of an external 
current at a given order in A but to all orders in as- This was done by including con¬ 
straints from collinear gauge invariance, the Dirac structure of the effective theory fields and 
reparameterization invariance. These conditions prove to be surprisingly predictive, and 
constrain not only the number of allowed operators, but also their functional dependence 
on label momenta. 

For the case of the heavy-light currents, the constraints from the Dirac structure of 
the effective theory fields have been included at leading order in and allow only 3 
structures in the = 0 frame. Here we consider the more general case of an arbitrary 
heavy quark velocity v, which is necessary in order to have a set of operators which closes 
under reparameterization transformations. 

At subleading order 0(A) the Dirac constraints alone allow many more operators. In 
particular, in addition to 2-body operators {^nW)iu ■ ■ ■ h^, one has to include also 3-body 
currents of the form (.CnhE)wi ■ ■ ■ (WHDW)^., ■ ■ ■ h„. RPI constraints on a subset of the two 


body operators were previously considered in 




and it was shown their coefficients are 


fixed in terms of the coefficients of leading order currents. Here we extended the constraints 
to the full set of allowed two-body and three-body operators, and showed that type (H) RPI 
imposed severe constraints on the (a;i,a; 2 ) dependence of the latter. For example, the scalar 
current qb is matched at 0(A) onto 8 general operators in the effective theory. After imposing 
all constraints, only one of these has a free Wilson coefficient, which has to be determined 
from a matching calculation. A similar reduction is obtained for the more complicated case 
of the vector/axial and tensor currents, for which one can write (28,44) structures but only 
(4,7) Wilson coefficients are not fixed by the symmetries of the effective theory. 

In this paper we have focused on mixed usoft-collinear interactions, however for many 
exclusive heavy-to-Iight processes the final operators that are needed are of soft-collinear 
type as was the case for heavy-to-light form factors. In practice it appears simplest to 
derive collinear-soft interactions from the collinear-usoft ones using the two-stage matching 
technique, QCD —> SCETj —SCETn, discussed in the proof of factorization for heavy-to- 
light decays Ref. Q. The operators in this paper describe interactions in the intermediate 
SCETi theory. For exclusive processes such as B ^ Dti where the intermediate ~ QA 
fluctuations in SCETj are responsible for inducing sim^e operators in SCETn the procedure 
used in Ref. Q reduces to the one discussed in Ref. HL 

Note Added: In the hnal stages of this work. Ref. [l2| appeared where a direct study of 
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light-light soft-collinear operators was performed. The Wilson coefficients of these opera¬ 
tors were determined by matching from QCD up to one-loop order, and both 2-body and 
3-body operators were found to contribute. An intermediate theory with modes of momen¬ 
tum ~ 1, A') as dynamical degrees of freedom was also considered. This appears 

similar in spirit to the QCD —> SCETi ^ SCETn construction used in Ref. 113 , however in 
the intermediate SCETj theory we found that the dynamical collinear modes should have 
momentum scaling as ~ Q(A^ 1, V^). Finally, reparameterization invariance constraints 
on soft-collinear operators were also discussed in Ref. and were shown to constrain the 
form of certain Wilson coefficients. 
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